Magneto-Resistance in the Two-Channel Anderson Lattice 
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The paramagnetic phase of the two channel Anderson Lattice model in the Kondo limit is investi- 
gated in infinite spatial dimensions using the non-crossing approximation. The resistivity exhibits 
a Kondo upturn with decreasing T, followed by a slow decrease to a finite value at T = 0. The 
decrease reflects lattice coherence effects in concert with particle-hole symmetry breaking. The 
magneto-resistance obeys an approximate scaling relation, decreasing towards coherent Fermi liquid 
behavior with increasing field. The magnetic field induces a Drude peak in the optical conductivity. 
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Introduction. Heavy fermion (HF) materials |l| have 
been under investigation for two decades. Character- 
istic for these materials is a 100 — 1000 fold enhanced 
electronic specific heat coefficient 7(T) = C(T)/T and 
very large and strongly temperature dependent resistivity 
p(T) (of the order of lOOpflcm). In most of these inter- 
metallic compounds magnetic or/and superconducting 
ground states are found. The physics of these interest- 
ing anomalous metals is related to strongly correlated 
electrons in 4//5/ orbitals. Some of these systems are 
well described as Landau Fermi liquids of massive quasi- 
particles. 

Recently, a number of non-Fermi liquid (NFL) HF al- 
loys have been found which display, e.g., logarithmically 
divergent j(T) 0. The superconducting HF compound 
UBei3 also has NFL behavior in 7, and moreover pos- 
sesses a very large residual resistivity (ss 100/ificm) at 
the superconducting transition even in high quality sam- 
ples (as determined by a large T c and sharp resistive 
transition). Such NFL behavior may be possible from 
a two-channel Kondo model description of the physics, 
such as has been proposed for UBei3 |3| on symmetry 
grounds. In this picture, electrical quadrupole moments 
of the two-fold non-magnetic T 3 ground state of the U - 
ion are screened by orbital motion of the conduction elec- 
trons. Because the magnetic moment of the electrons is 
a spectator to this process, there are two screening chan- 
nels. Reversal of spin and orbital indices allows for a 
two-channel magnetic Kondo effect for a Ce 3+ ion in a 
cubic environment [Q . The two channel Kondo impurity 
model for these cases |||| has been investigated essen- 
tially exactly with different techniques However, lit- 
tle is known about the corresponding lattice model . 

In this paper we present a solution of one particle prop- 
erties of the two channel Anderson lattice model in in- 
finite spatial dimensions. To obtain the solution of the 
effective two channel single impurity problem, the Non- 
Crossing-Approximation (NCA) ||-|ll],D is use d- Al- 
though this method fails to solve the infinite dimension 
single channel Anderson lattice at temperatures less than 
the lattice Kondo-temperature T* jl2|, it works well in 



the low temperature region of the two-channel model 
fjl,[i3p|l . The calculated resistivity p(T) agrees well with 
recent Quantum- Monte-Carlo (QMC) results ||. For 
lower temperatures than accessible by QMC, however, 
we find a decrease in p(T) for decreasing T. A strong 
negative magneto-resistivity at temperatures below the 
resistivity maximum in an applied magnetic field indi- 
cates a recovery of Fermi-liquid behaviour. At the same 
time, an onset of a Drude peak in the optical conductiv- 
ity is found which is absent in the zero field solution. We 
briefly discuss the possible relevance of these results to 
experiment at the end. 

Theory. The two channel Anderson lattice Hamilto- 
nian under investigation reads 
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X are the usual Hubbard operators, d being the spa- 
tial dimension, i the lattice site, t* the reduced hop- 
ping matrix element of the conduction electron between 
nearest neighbours which carry a spin a and a channel 
index a = (1,2). The conduction electrons couple via 
hybridization matrix element V to the ionic many-body 
states on each lattice site. The symmetry breaking mag- 
netic field enters by a Zeeman term in E a — Sf+g\pb\o"H. 
The Zeeman splitting of the conduction electrons only re- 
sult in a shift of the band centers and turn out to be a 
small correction. When \E a — E a \ is much larger than 
the hybridization width Tq = irpV 2 , the model can be 
mapped onto an two channel Kondo model |3| via the 
Schrieffer- Wolff transformation |Q . 

In the local approximation for the Anderson Lat- 
tice [Q, which is equivalent to the limit d — > 00 flq ] 
with an appropriate rescaling of the effective hopping 
]l6[ , we choose one /-site as an effective '"impurity"' 
site which is self-consistently embedded in an effective 
medium reflection the rest of lattice. While in a sin- 
gle impurity problem only the bare medium A^ a (z) = 
N ^ Efc d aa {k, z) = N- 1 J2 k y2 ( z ~ £fea<r) _1 enters, in 
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the lattice the condition that the local /-Green's function 
(GF) has to be equal to the /c-summed lattice GF 
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where A aa (z) is the self-consistent one body self en- 
ergy of the local "impurity" propagator F aa {z) =<C 
Xa]a\Xa}a 3> (z). The effective hybridization width 
r a(T (w) = ^K atT {ui — iS) enters then the local two chan- 
nel impurity problem. The self-energy of the conduction 
electrons is given by 



V 2 F a(7 {z) 
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In the single channel case, the unitarity limit of the 
T-matrix T(z) = V 2 F(z) leads to the Fermi-liquid be- 
haviour of the conduction band self-energy T, a (z) oc 
aT 2 + bw 2 . Since the value of the T-matrix at the chem- 
ical potential and T — > is smaller than the unitarity 
limit in the two(multi)-channel case Eqn.(§) tells us im- 
mediately, that the corresponding conduction band self- 
energy for the exact solution in the paramagnetic phase 
of the lattice has to be finite. This has been recently 
called an incoherent metal Q]^]. The physical origin is 
the following: the local spin is over-compensated by two 
conduction electron spins. On each lattice site a resid- 
ual free thermodynamically fluctuating degree of freedom 
(DOF) acts as a scatterer for conduction electrons. An 
residual entropy of l/21og(2) per site is associated with 
this DOF freedom which has been interpreted as a free 
Majorana fermion Jx7[ ■ The finite self-energy yields a fi- 
nite value for p(T — > 0, H = 0). Since in translational in- 
variant system a vanishing dc resistivity is expected this 
indicates that the paramagnetic state is not the ground 
state of the two-channel Anderson lattice. 

In the absence of a magnetic field the NCA equations 
of the effective impurity 
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are equivalent to a resonant level system with an effective 
Anderson width f = 2r . The NCA pathology in the 
local GF becomes the physical Abrikosov-Suhl-resonancc 
(ASR) in the two channel case The NCA threshold 

exponents of the effective ionic propagators P |lq] have 
the exact value of 1/2 calculated within a conformal field 
theory approach jl9| |. In limit of infinite spin N and 
channel M degeneracy with a fix ratio N/M the NCA 
becomes exact |l3j. The effective local GF is given by 
the convolution 



1 f dz 

F aa {iu n ) = -~-i —e- pz P a {z)P a {z + iw n ), (6) 
Zf Jc 2m 

where Zf is the effective local partition function. Even 
though higher order vertex corrections pC( ] will modify 
the spectral distribution, the leading physical effect and 
the correct thermodynamics is captured correctly within 
the Eqns.(^) and . Since the saturation value of the 
effective site T-matrix is half the unitary limit no pseudo- 
gap develops in the quasi-particle spectrum as in the one- 
channel lattice jH) . 

In infinite spatial dimensions the vertex corrections in 
the two-particle propagators vanish [^2| and the conduc- 
tivity itself is a 1/d correction which can be calculated 
by evaluating the lowest order bubble diagram f23| , given 
by 

OO oo 
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which can be written as an integral over four complex er- 
ror function; A = ne 2 a 2 t* 2 N(hdVol)- 1 = t* u? p / '(4tt) , the 

Gaussian density of states po(e) |l6|, Gal{z) the con- 
duction electron GF, and a the lattice constant of the 
d-dimensional hypercube. The /-electrons do not con- 
tribute to the conductivity since the hybridization is k in- 
dependent. The dc-conductivity is obtained by the limit 
o-dc(T) = lim a(uj,T). 

Results. To obtain a self-consistent solution of the lat- 
tice problem, (i) the effective hybridization with r(e) has 
been treated with the same accuracy as the threshold sin- 
gularity of the ionic propagators in Eqns.(^) and (||), and 
(ii) only 10% of the calculated change of r(e) is added 
in each lattice iteration step, i.e. the lattice is switched 
on adiabatically. The error in norm of the P(z) reaches 
0.01%, the sum rule for the self-energy is obeyed within 
0.02% and the maximum norm max{\T n (e)— r n _i(e)|} < 
10~ 8 . All energies, if not otherwise stated, are mea- 
sured in the original Anderson- width Tq. We chose 
Ef = E a — E a — — 3r in the absence of H and t* = 10r 
with a band center at uj = 0. T K = 0.016r |||. 

In Fig. HI p(T, H)) normalized to the estimated T — > 
value of the QMC data || is shown for different val- 
ues of the an applied magnetic field measured in units 
of Hk = ksTx I (<?Mb)' We have fix the lattice scale 
T = 1.3Tk by matching the QMC resistivity data. 
The agreement with the higher temperature data for the 
Kondo lattice - the open symbols - is excellent. Nev- 
ertheless, the resistivity has a maximum and slowly de- 
creases with decreasing temperature, as expected from 
a lattice calculation. The calculations are done slightly 
below half-filling and with fixed chemical potential. At 
half-filling in the Kondo-regime an the analytical solution 
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obtained with an artificial Lorentzian density of states |t]] 
predicts a resistivity of p(T) = p(0)(l — aVT) at low tem- 
perature. While for H = still a positive intercept at 
T = is expected consistent with an infinitly degener- 
ated ground state, clearly a crossover to Fermi-liquid is 
found within an applied field. The NCA pathology Q 
prevents access to the Fermi-liquid regime when T — ► 
pjfl . Evaluation the constants for A in Eqn. (J?]), assum- 
ing a lattice constant of 5A and 2 electrons per unit cell 
in a three dimensional lattice gives a resistivity prefactor 
of « 12.6fiflcm/TQ which leads to a resistivity maximum 
of w 250/i£lcm using our absolute maximum of 20Fq. 
This is very close to the experimentally found value of 
« for UBei 3 ||. 

Motivated by the experimental data for p(T, i?) for 
UBei3 p7f| , we have attempted to scale our p(T,H) data 
with the Unsafe Ap/p = [p(T,H) - p(T, 0)]/p(T, 0) cx 
f(H/(T + T*Y). While for the impurity model, we ex- 
pect T* = 0.,/3 = 1/2, we find approximate scaling for 
T* = 0.006T^,/3 = 0.39, as ploted in Fig. |. The inset 
of the figure shows the imaginary part of the conduction 
electron self-energy E c (w — id) (§) for H = is plotted 
for four different temperatures in Fig. 0. It shows a shift 
of the maxima away from the chemical potential in this 
metallic regime. Very close to uj — a very small onset 
of coherence is observed for T — > 0, but the relaxation 
rate remains of the order of 2ro 

Now we focus on the optical conductivity, displayed 
in figure m The large peak at w 0.9r results from high 
energy charge excitations. With decreasing temperatures 
the optical conductivity develops a pseudo-gap. The /- 
sum rule relates the integrated optical conductivity 
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to the average kinetic energy in the direction of the cur- 



rent flow 1 26 1 , which is checked numerically. It indicates a 
shift of small amount of spectral weight to higher frequen- 
cies. This can be seen clearly in the figure by comparing 
the T = 10T K and the T = T K curve (note that the log- 
arithmic plot overemphazies the area of the gap) . At low 
temperatures a small increase in cr(u>) can be observed 
when u> — > 0. Nevertheless, no clear Drude peak is seen 
even for T = O.OlTff, one decade lower than the observed 
maximum in 1/<t(0, T). However, in a magnetic field of 
H = Hx a low frequency "Drude" -peak develops again, 
consistent with the return to Fermi-liquid behaviour sug- 
gested in p(T, H). Note that in the single channel Ander- 
son lattice a clear Drude peak develops already a little 
below the maximum in l/tr(0,T). 

Using the Kramers-Kronig relation the imaginary part 
of the optical conductivity has been calculated. With the 
phenomenological Ansatz 
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the dynamical optical relaxation rate Y opt (<jj) and the 
dimensionless mass enhancement factor A(u>) have been 
determined. In Fig. ||(a) r opt (cj) is plotted for the 
same parameters as in Fig. |3|. While for temperatures 
T > OATk T opt {uj) is nearly frequency independent for 
low frequencies, at the lowest temperature r opt (0) has 
decreased reflecting the decrease of the dc-resistivity. 
The low frequency behaviour has an exponent slightly 
lower then n = 1. We emphasize the difference between 
3r7iE c (Li;) and T opt (w): the first is a true one-particle 
relaxation rate, the second, however, reflects the two- 
particle nature of the energy absorption process associ- 
ated with electrical charge transport. Generally, only for 
a Fermi-liquid at very low temperatures and frequencies 
should T opt = -23mS c (w). In an applied magnetic field 
H = Hk and low temperatures T — O.OITr- the opti- 
cal relaxation rate r opt (0) shows the expected trend to 
Fermi-liquid behaviour. 

Fig. ||(b) shows the quantum Monte Carlo (QMC) cal- 
culation of T op t(u>) for the two-channel Kondo lattice at 
particle hole symmetry. A strict quantitative calculation 
is not possible because: (i) of the overlap of To and the 
Kondo interaction J to within an order of magnitude in 
the QMC calculations (in the Anderson model calcula- 
tions of this paper Tq is well separated from high energy 
scales), and (ii) because the particle-hole symmetry re- 
moves the non-monotonicity experienced in the NCA cal- 
culations. Modulo these concerns, the separate calcula- 
tions agree qualitatively in the overlapping temperature 
and frequency regions. 

Comparison to Experiment: As mentioned earlier, 
p(T, H = 0) is reminiscent in form and magnitude to 
UBei3 P9,E7|. The magneto-resistance also resembles 
that of UBei3, though our scaling form in detail is differ- 
ent. However, a strict comparison is not possible, since 
assuming a quadrupolar Kondo model applies to UBei3, 
we should rather split \ai > states (order H) and quadrat- 
ically split \ai > stated (van Vleck processes). Addi- 
tionally, recent experiments suggest a possible U 3+ -U 4+ 
configuration degeneracy with is lifted with Th substi- 
tution p8| . We refer the necessary intermediate valance 
calculation to a future work. In this case, a crossover 
from NFL to Fermi-liquid physics is still expected. Our 
cr(T, ui), T opt {ijj) calculations are very compatible with 
data for the alloys Y0.sU0.2Pd3 and Thi^U^Pds Al 3 [§| , 
as well as the compound UBei3 (3^] . Because of the in- 
coherent normal metal phase, we expect little qualitative 
difference between these more dilute alloys and the lat- 
tice. For UBei3, the existing optical data only go to 50 
cm -1 ~ 5 — &1zbTk in frequency p0| , and it is clearly de- 
sireable to extend these measurements to lower frequen- 
cies. We remark that for Thi_ a; U 2 ;Pd2Al3, if a hexago- 
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nal quadrupolar Kondo picture applies, a c-axis magnetic 
field will split the \ai > levels S, permitting comparison 
to our calculations. In addition, it should be interesting 
to test whether a magnetic two-channel lattice picture 
applies to CeCu 2 Si2 [Q and thus have detailed optical 
conductivity measurements carried out in applied field 
for this system. 
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FIG. 1. Resistivity for the two channel Anderson lattice model (TCA) vs. temperature for different magnetic fields. We have 
normalized to the estimated p(T = 0) values of the Quantum-Monte-Carlo (QMC) data, temperature for the same parameters. 
The open symbols are the QMC results for different J in the two-channel Kondo lattice model. The inset shows the scaling of 
the position of the dc maxima vs H. Parameters: e/ = — 3ro,i* = I0ro,Tb = 0.02ro. 
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FIG. 2. Magneto-resistance {p{T,H) - p(T, 0))/p(T, 0) vs scaling variable x = H/(T + 0.006Tk ) ' 39 . Inset: Imaginary part 
of the conduction band self-energy vs. frequency in the vicinity of the chemical potential /i = for four different temperatures. 
Parameters: ej = —3To,t* = IOFq. 
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FIG. 3. Optical conductivity per channel in units of Wp/(47r) vs. frequency for covering three decades in temperature. In 
the inset the corresponding curves for the single channel PAM for the two highest temperatures are show in the same units. 
The dash-dotted line is calculated with H — Hk- Parameters: £/ = — 3Po,t* = IOPq. 
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Two Channel Anderson Lattice 
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FIG. 4. Relaxation rate F op t(u}) in units of To (a) per channel for the two channel Anderson lattice measured versus lj/Tk 
(Parameters: £/ = — 3To, t* = IOFq) and (b) QMC data for the two channel Kondo lattice (Parameter: J = 4Fo, To = 0.281Fo). 
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